Metastability and uniqueness of vortex states at depinning 
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We present results from numerical simulations of transport of vortices in the zero-field cooled 
(ZFC) and the field-cooled (FC) state of a type-II superconductor. In the absence of an applied 
current 7, we find that the FC state has a lower defect density than the ZFC state, and is stable 
against thermal cycling. On the other hand, by cycling I, surprisingly we find that the ZFC state is 
the stable state. The FC state is metastable as manifested by increasing / to the depinning current 
7c, in which case the FC state evolves into the ZFC state. We also find that all configurations 
acquire a unique defect density at the depinning transition independent of the history of the initial 
states. 

PACS numbers; 74.60.Ge, 74.60.Jg 



Metastability is a generic feature in many disordered 
systems. Metastable states show characteristics distinct 
from the ground state, and are evidenced, for example, 
by the path or history dependence of the different re- 
sponse functions. A striking example of metastability 
is the difference between the field-cooled (FC) and zero- 
field cooled (ZFC) magnetisation in vortex phase of type- 
II superconductors Ij. In superconductors, the FC mag- 
netisation is much smaller than the ZFC magnetisation. 
This difference decreases with increasing temperature T 
and vanishes above the irreversibility temperature Tirr- 
The ZFC state is believed to be metastable and the FC 
state to be a quasi-equilibrium state because it is possible 
to reach the FC state from the ZFC state by changing 
T from a low value to Tirr and back, but the ZFC state 
cannot be reached from a FC state through such a tem- 
perature cycle. Similar features are also observed in spin 
glasses The difference in FC and ZFC state magneti- 
sation is explained by the emergence of energy barriers 
below Tirr, which inhibits the ZFC state from exploring 
the full configuration space. 

In recent years, transport experiements have shown 
conflicting results regarding the stability of the FC and 
ZFC states. In the anomalous peak effect region, the 
FC state has been observed to be unstable against cur- 
rent annealing. The critical current Ic is higher (upto 
a factor of 6) in the virgin FC state than in FC states 
annealed by a current or in the ZFC state 0]. On the 
other hand, using a slow ramp rate for /, the FC and 
ZFC state Ic was found to be same below and inside the 
peak effect regime^. This was interepreted as a sig- 
nature of the metastability of the FC state. Also, the 
susceptibility of the FC state can be switched to that 
of the ZFC state by applying an ac pulse 0, • More- 
over, recent scanning Hall probe measurements have also 
shown substantial vortex rearragement in FC states for 
I < Ic 0- One approach accounts for these experimen- 
tal observations through the role of surface barriers 
These differing pictures place in the focus the behavior 
of vortices in the FC and ZFC states in the presence of 



a transport current. 

In this paper we use numerical simulations to study 
the transport of differently prepared vortex systems. Our 
first main result is that, in contrast to the temperature 
cycling, if the vortex system is cycled with a transport 
current, then the ZFC state is more stable than the FC 
state. In particular, the ZFC state can be reached from 
the FC state by current cycling, but not vice versa. Our 
second main result is that, with increasing currents, all 
vortex states, FC, ZFC, or prepared in any intermediate 
way, acquire an unique density of defects at depinning. 

We consider vortices in a 2D plane perpendicular to 
the magnetic field B = Bi, of a bulk superconductor, 
with B = ny(j)Q where is the vortex density, and (j)o is 
the flux quantum. Within the London approximation we 
treat the vortices as particles, with dynamics governed 
by an overdamped equation of motion, 

j^i k 

Here rj is the flux-flow viscosity, and the first term rep- 
resents the inter-vortex interaction given by the poten- 

tial U'"{r) — g^i^-i Ko{r/X), where Kq is the zeroth- 
order Bessel function, and r = {r'^ + 2^^)^/^ with A 
and ^ as penetration depth and coherence length of the 
superconductor, respectively 0- The second term is 
the attractive interaction with parabolic potential wells 
[/p(r) = Uo{^ — 1) for r < Vp, and otherwise, cen- 
tered at the random locations. The third term 
Fext = ■jJ X (j)QZ represents the Lorentz force due to the 
transport current density J. The last term is the thermal 
noise (Langevin term) with the moments (Ci,p(i)) = 0; ^^nd 
{Ct.p(t)Q,p'{t')) = 2kBTr]S^jdpp'S{t - t'), where T is the 
temperature, ks is Boltzmann constant, and p,p' = x,y. 
We worked with reduced variables: all distances are in 
units of Aq = A(r = 0), the current density J in units 

of c/o/00, and T in units of Xofo/ks where /o = g^^- 
The time t is in units of ?/Ao//o. This model is expected 
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to reasonably capture the physics of stiff vortex Hnes, 
hke those in Nb and NbSe2, while in strongly layered 
superconductors, the soft interlayer couplings introduce 
additional degrees of freedom. 

To simulate the ZFC state, we confine the quenched 
disorder to a central region of the simulation box, while 
leaving the outer region defect free as illustrated in 
Fig. 1. The disordered region defines the extent of the 
superconductor and the disorder free region simulates 
the free-space. Periodic boundary conditions connect 
the edges of the simulation box. We employ a box of 
size 40 Ao x 24 Aq, with the disordered region occupying 
28Ao X 24Ao. We used a pin density of 5.95/Aq with 
(t/o) =0.03, and = O.IAq. The results presented below 
are for Ny = 1148 vortices in the simulation box which 
corresponds to B fa 800G. Typically, we used 5 x 10"* 
time steps for thermal equilibriation. For V-I character- 
istics, 2 X 10"* time steps was used to equilibriate before 
averaging over a similar time scale. The real space config- 
uration of vortices is characterized by the defect density 
Ud — Nd/S where Nd is the number of voritces which 
are non six-fold co-ordinated and S is the area (in units 
of Aq) of the disordered region 0. A^^ is counted by 
Delaunay triangulation of the real space position of vor- 
tices. The V-I characteristics are calculated by varying 
^ext = Fext^ {Fext oc the Current) and calculating the 
average velocity (oc V, the voltage) v = Sf" ^^ong 
the a;-axis. 

We first consider the effect of cycling the temperature 
T. The ZFC state at T — is created by positioning vor- 
tices in the defect free region and allowing them to move 
into the disordered region by their own dynamics. The 
FC state is created by the thermal cycle of the ZFC state 
by first heating the system deep into the liquid phase, 
and subsequently cooHng it back to T = 0. Fig. 1 shows 
the behaviour of the defect density Ud throughout this 
thermal cycle. The ZFC state at T = is formed by vor- 
tices diffusing through the free edges into the disordered 
region. With decreasing velocity of the diffusion front, 
plastic deformation induces a larger density of defects. 
These defects are non-equilibrium and is supported by a 
gradient (oc (t/^)) in vortex densitv[ll|. On increasing T, 
the gradient, and subsequently n^(T), decreases due to 
thermal activation. Across T = 0.012 ^T^, the melting 
temperature, nd{T) increases rapidly as the system goes 
into the liquid phase. Compared to the ZFC state, the 
FC state has a lower nd{T — Q) as it is cooled from the ho- 
mogeneous liquid phase that has no macrosopic density 
gradients. The FC state shows no appreciable change in 
Ud on repeated thermal cycling. The inset in Fig. 1 shows 
the hysteresis in Ud on crossing Tm, which is evidence for 
the melting transition into the liquid phase being first 
order. A very slow cooling rate is used in the inset com- 
pared to the main panel to more accurately observe the 
hysteresis across T^. 

A partial thermal cycling of the system, consisting of 




FIG. 1: Plot of nd{T) for the thermal cycling starting with 
the ZFC state at T = 0. Also shown is the partial stability of 
the PFC state, created by cooling from Tpfc = 0.006. The 
direction in which T is changed for creating the PFC state is 
shown by thick arrows. The simulation geometry is shown on 
the upper left corner of the main panel. Inset : hysteresis in 
Hd across Tm during the FC cooling and warming procedure. 

a warming step to a temperature Tppc, e.g. Tppc — 
0.006 = 2iL^ and a subsequent cooling step to r = gen- 
erates a partially field-cooled (PFC) state (see Fig. 1). 
The PFC state is less disordered than the ZFC state, 
and is stable to thermal cycling between and Tppc- 
For T > Tppc, the PFC state follows the same path as 
the virgin ZFC state. This repeated warming and cool- 
ing allows one to access all the other PFC configurations 
at r = 0. The behaviour of nd{T) is in consonance with 
the experimental measurement of magnetisation M{T) 
for T<T,rM. 

Next we consider the effect of current / on various con- 
figurations. As shown in Ref.|l3|, effects of the disorder 
on the moving vortex state can be characterized by a 
"shaking temperature" Tsu- Tgh oc v~^, where u oc / is 
the average vortex velocity. At large u, a moving lattice 
is formed which undergoes an order to disorder transition 
as / is decreased to zero across the depinning current Ic- 
Therefore, besides the above described thermal cycling, 
a disordered vortex state can be prepared by current cy- 
cling as well. As will be shown below, the resulting con- 
figuration is the most stable vortex state against repeated 
cycling of /. We call the resulting static configuration 
the current hardened (CH) state, in analogy to "strain 
hardening" , in which materials are cycled with increasing 
stress to increase the density of dislocations. 

Fig. 2 shows V{I) and nd{I) curves at T = for differ- 
ently prepared states. We identify three current regimes 
: (1) a pinned regime I<Ic, with V = Q and nd{I) deter- 
mined by the mode of preparation; (2) a plastic regime 
Ic<I<Im, with non-linear V{I) and the same nd{I) for 
all differently prepared states; and (3) a flowing-lattice 
regime with linear V{I) and nd{I) substantially 
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smaller. The small difference in Ud in the flowing-lattice 
regime, as seen in Fig. 2, is due to the formation of a grain 
boundary at the free edges of the disorder region (see also 
Fig. 4(c)). Remarkably, V{I) curves and the depinning 
current Ic are independent of the mode of preparation of 
the system in agreement with the experiments with 
slow ramping of These characteristics persist at fi- 
nite T as well. The nd{I) for I < Ic, on the other hand, 
depends on the mode of preparation. 




FIG. 2: V{I) (filled symbols), and nd{I) (open symbols) at 
r = for the FC state (o), PFC state (v), ZFC state (□), 
and CH state (A). The direction in which current is changed 
is shown by thick arrows. Inset : same as the main panel for 
r = 0.004. The dashed and the dotted lines mark lirr and Ic, 
respectively. 

Strikingly, the defect density nd{I) at the depinning 
current Ic assumes the same value n^sa 0.51 for the FC, 
ZFC, PFC, and CH states, as seen in Fig. 2. In other 
words, is independent of the mode of the preparation 
of the vortex state. The history independent unique value 
of implies that the characteristics and dynamics of all 
vortex states converge as the current approaches Ic- This 
idea is substantiated by studying the real space dynam- 
ics at depinning. As is known, vortices depin by forming 
few channels across the system [3- We find that, while 
in differently prepared states the channels are nucleated 
at different locations along the entry edge, in the bulk 
the channels converge quickly to a universal pattern. We 
also find the average transverse wandering of a channel 
is of the order of ~ 6-7 lattice constants, and is same for 
all states at Ic- This implies that the average displace- 
ment of the vortices at Ic is the same for all states, thus 
providing an underlying physical picture for preparation 
independence of n^. Further studies show that the value 
of nj^ depends on T, n^, and the disorder strength. 

Next, we study the stability of the different vortex 
states against cycling with I. The virgin FC state has 
a defect density nd,Fc{0) <n'^- Therefore, the density of 
defects in the virgin FC state increases from nd,Fc(0) to 
n"^ as the / is increased from to Ic, as shown in Fig. 2. 



Furthermore, decreasing current from Ii < Ic, n-d.Fcil) 
does not decrease to its virgin value, but stays approxi- 
mately at nrf(/i). In particular, cycling / between and 
I>Ic leaves nd,Fc{I) ~ "-^ ^'^^ I<Ic- At the same time, 
nd,ZFc{0) ~ n'^, and cycling / does not affect this value. 
Therefore, the FC state, and its defect density Ud.FC, 
can be transformed into the ZFC state and Ud^zFC by 
cycling / between and Ic- The reverse is not true, i.e., 
the FC state cannot he reached from any state by cycling 
current- Therefore, in contrast to thermal cycling, we 
find that the FC state is metastable and the ZFC state is 
stable against current cycling. Fig. 2 further illustrates 
that the CH state and its defect density nd.cH{I) evolve 
into the ZFC state when the current is decreased from a 
large value to 0. The same holds for the PFC states, as 
nd,PFc{I) also evolves into Ud.zFci^)- 

Since for I > Ic nd{I) is approximately the same for 
all states, metastable configurations and the irreversible 
behavior of nd{I) emerge only below Urr- For T — Q, 
lirr = Ic, but at finite temperatures these two current 
values separate, as evident for T — 0.004, shown in the 
inset of Fig. 2. At finite T, nd{I) « nj^ at a current Urr > 
Ic- This is possibly due to thermally activated vortex 
creep, which induces vortex displacements already for / < 
lirr- The metastability of the FC state persists at finite 
T, and we believe is a generic phenomenon which brings 
out a fundamental difference in the nature of FC and ZFC 
states. The ZFC and CH states show similar behaviour 
throughout the studied current and temperature ranges. 
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FIG. 3: (a) The time evolution of rid at T = 0.0 for a pulse 
of current I = 0.15 < Ic in the ZFC and the FC state, (b) 
Plot of Ic{T) obtained from FC, ZFC and CH states, (c) 
Shows n'^{T) obtained from V{I) curves along with Ud for 
ZFC warming and FC cooling superimposed from Fig. 1. 

Further evidence for metastability of the FC state is 
obtained by applying a current pulse I — 0-15 < Ic to dif- 
ferently prepared vortex states, and observing the tem- 
poral evolution of Ud, as illustrated in Fig. 3(a). The 
application of a current pulse evolves the FC state into 



PFC states with an intermediate defect density nd,FC < 
nd,PFC < nd,ZFC- The defect density remains at that 
elevated level even after the current drops to zero. In 
particular, a large enough current pulse evolves the FC 
state into the ZFC state. This transformation of the FC 
state to the ZFC state by a currentpulse is in close anal- 
ogy to experimental observations [3 ■ In Ref. 0| , the Ic is 
also found to change during the transformation of the FC 
state to ZFC state, possibly due to the influence of sur- 
face barrier TBI . We emphasize that the transformation of 
the FC configuration to ZFC configuration as found here 
is a bulk characteristic, and is not influenced by the free 
edges of the disordered re gion . This was verified in sim- 
ulation without free edges The FC state was again 
found to be metastable below /irr, and is transformed 
into CH configuration on increasing / above hrr- 

In Fig. 3(b) we plot /c(T) for a range of T at which 
simulations were carried out. The Ic follows an exp(— ^) 
form, confirming that thermally activated motion domi- 
nates the dynamics at Ic- Fig. 3(c) shows the tempera- 
ture dependence of n'^{T), as obtained from V{I) curves. 
The nd{T) values for the FC and ZFC states are super- 
imposed from Fig. 1. n^{T) closely follows nd,zFc{T), 
obtained by warming the ZFC state, thus corroborat- 
ing the view that the ZFC state is essentially the frozen 
equivalent of the critical state at Ic- In other words, the 
channels formed by the penetrating front of vortices dur- 
ing the formation of the virgn ZFC state is frozen. Any 
subsequent application of / essentially uses these pre- 
formed channels which explains why ZFC state shows no 
substantial change in Ud- In contrast, in the virgin FC 
state no such channels exists. On first passage of /, the 
resulting build up of stress at depinning leads to chan- 
nel formation which increases Ud- Any further cycling of 
/ essentially uses the already formed channels. This is 
a plausible explanation for the metastability of the FC 
state on cycling /. 

The configurational changes are illustrated in Fig. 4 for 
various states, during a composite T-I cycle. Starting 
with the ZFC state, shown in Fig. 4(a), we obtain the 
FC state (Fig. 4(b)) by cycling T from to T>T„. The 
FC state is then driven with / > /,„ to obtain Fig. 4(c). 
And finally, decreasing / to 0, the CH state is obtained in 
Fig. 4(d), which is indistinguishable from the ZFC state 
in Fig. 4(a). Thus a composite T ~ I cycle allows one to 
go from the ZFC state to the FC state and vice-versa- 

In conclusion, we have shown that the ZFC state is 
stable against current cycling in contrast to the FC state 
which is found to be metastable. This behavior is oppo- 
site to the thermal cycling where the FC state is stable. 
We also find that the defect density Ud for any initial 
state converges to the same value at Ic- 
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FIG. 4: The real space configuration at r = during a com- 
posite T-I cycle : (a) ZFC state, (b) FC state, (c) flowing 
lattice state, (d) CH state. The arrows show the cycling path 
(see text for details). Note the close similarity between (a) 
and (d). The symbols (o) and (•) denote vortices with 5 and 
7 neighbours, respectively. 
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